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eit CDUCTION 


An overall rise in fuel prices has led to an increasing 
interest in the design of autopilots for ships. The purpose 
of the automatic steering control is to minimize propulsion 
losses which are caused by added drag due to steering of the 
ship. 

Several attempts have been made to define a measure of 
steering efficiency based on propulsion losses and proposals 
have been made for the design of an autopilot which mini- 
mizes the propulsion losses. As a performance criteria, 
added resistance due to steering suggested by Norrbin, N. H. 
was used in most of the studies because it is convenient for 
ship board use. 

Many researchers indicate that a carefully designed 
controller could save from one to two percent of fuel by 
minimizing the propulsion losses which are caused by added 
drag due to steering of the ship . No analytical evaluation 
of propulsion losses due to steering ina sea state has been 
made. 

The goal of this thesis is to study fuel savings using 
various kinds of controllers, and especially to compare an 
Adaptive Automatic Optimal Controller with a PID controller 
to see the difference between them. The optimal gain parame- 
ters of the Adaptive Automatic Optimal Controller are 
provided by LTJG. Cetin Diken [Ref. Ale 

To study the optimization problem, models of both the 
ship and the operating environment are required. Chapter 2 
addresses what type of computer model can be used to repre- 
Sent the ship. 

The ship's nonlinear equations of motion were needed to 
Simulate the ship in the computer program. Chapter 3 
addresses the Mariner Class ship nonlinear equations of 


motion. 


Bat 


Chapter 4 addresses the problems of testing the ship 
Simulation model in calm water, of expressing the forces 
and moments due to sea state, and the effects of sea state 
on ship's behaviour. 

Chapter 5 addresses the derivation of the cost function 
which represents the added drag due to steering. 

What is an adequate way to represent the fuel consump- 
tion ratio of added drag due to steering? Chapter 6 
addresses the fuel consumption ratio. 

Ship dynamics change with operating conditions such as 
ship speed, encounter angle, encounter frequency, and sea 
state. Chapter 7 presents the fuel consumption ratio of 
added drag due to steering for PID and an Automatic Optimal 
Controller, and the difference between them. 

Conclusions were drawn from simulation results, and are 
presented in Chapter 8. This chapter also recommends topics 


for future studies. 
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ize PeecCktenom OF Tht COMPUTE ReMODEL 


Before a specific controller can be designed, a real- 
Hocic model for the ship dynamics must be found to enable 
proper simulation of the ship movements in response to 
Bomcrol signals. 

The model which best represents ship-steering dynamics 
is a Taylor's series expansion of the force and moment rela- 
tionship around a selected steady-state operating point. The 
resulting equations are commonly Known as the equations of 
motion.{Refs. 2, 3, 4, 5] 

A computer program was developed in order to provide a 
computer simulation for the ship, using experimentally meas- 
ured hydrodynamic coefficients for the Mariner Class Ship 
(Refs. 2, 3, 6, 7}. Figure 2.1 shows the block diagram. The 
Bemeuter program is shown in Appendix C [Refs. 8, 9, 10, 
ie 

The function minimization subroutine used was BOXPLX 
which was programmed by R. Hilleary. The task of BOXPLX is 
to find the minimum of any function. It may be subjected to 
meoricit Constraints of the variable or implicit constraints 
on functions of the variables. [Refs. 11] 

The sea disturbance is found as first and second-order 
forces and moments by a sea state program which has been 
written by J. Cass [Ref. 9]. 

D type or PID controller was used in the computer simu- 
lation program. PID controlller is shown in Appendix A. The 


D type controller was described by: 


GI Ibs Hee i 


(1+T,s) ee 


3 


where I le ao ame! aS are optimal parameters which were 


One 
found by the minimization subroutine, BOKPLXA. The cena 
tions of symbols used in equations and figures are given in 


Appendix D. 





Figurewz a. Block Diagram of Ship and Contra 22 


ITI. NONLINEAR EQUATIONS OF MOTION 


Linear theory is useful for analyZing the influence of 
ship features on controls-fixed stability as well as on the 
turning ability of stable ships in the linear range [Refs. 
Pee ao, 4, Se l2Z2, 13). However, it fails to predict accurately 
the characteristics of the tight maneuvers that most ships 
are capable of performing and it cannot predict the maneu- 
vers of unstable ships. 

Nonlinear equations ous motion are suitable i @ 3 
predicting tight maneuvers and also suitable for computer 
programming. The nonlinear equations of motion based on abe. 
“Taylor series expansion of forces and moments including 
terms up to the third order have been developed by AbkKowitz 
[Refs. 5, 12], and Strom-Tejsen [Ref. 3]. Accuracy is not 
improved by including terms higher than third order. The 
development of these nonlinear equations is based on a 
restatement of linear equations to include rudder angle. 
Peiactons X, Y, and N are functions of u, v, Y, U, V, F, and 
Q - 

Combining the nonlinear Taylor series expansion of 
forces and moments terms up to third order with the dynamic 
response terms of the X, Y, and N equations, the nonlinear 
equations of motion are shown in Appendix B. The hydrody- 
namic coefficients of the equations are determined from 
experimental data obtained from captive model tests.They are 
given in Appendix B [Ref. 7]. The terms not included in 


tables are negligible. 


> 


IV. SHIP'S BEHAVIOUR IN CALM WATER AND SEA STATES 


A. CALM WATER 

A simulation program in Appendix C was run for turning 
and zlg-zag maneuvers to observe u, v, fr, and w using the 
ship's nonlinear equations of motion and Mariner Class Ship 
coefficients. 

It was observed that the rudder angle changes the ship's 
course. The ship's speed decreased while turning. The abso- 
lute value of v and r increased and after ‘a while they 
reached steady-state values, since the ordered rudder angle 
is constant. The larger the rudder angle, the greater the 
decrease in speed, and the faster the ship turns. 

Time responses of r, v, and u are shown in Figure 4.1, 
4.2, 4.3, and X vs. Y is shown in Figure 4.4 for the turning 
maneuver with 9 = 25 degree. §, Ww, and u are shown ini age. 
4.5 and Figure 4.6 for the Zig-Zag maneuver. As is seen from 
the figures, a suitable and sufficiently accurate ship 


computer model was defined. 


B. sbAgeodaik 

To observe the ship's behavior in a sea state, distur- 
bance forces and moments are needed. They depend on sea 
state, ship speed, encounter angle, and encounter frequency. 
The added mass and added inertia are functions of encounter 
frequency and sea state. 

A regular sea model was used as the sea representation: 
the wave crests assumed to be straight, infinitely long, 
parallel, and equally spaced with constant wave height 
[Refs. 2, So ia yeeros 

The forces exerted by the regular sea can be represented 
as the sum of two components, called first and second-order 
forces [Refs. 2, 6, 9, 14, 16, 17 )138)3c oe 


iG 


The first-order exciting forces are a zero mean process. 
The second-order forees are a nonzero mean process, il1.e. , 
they have essentially constant values. In regard to maneu- 
vering, the first-order forces are not the primary mover of 
the ship. The reason for this is that the high frequencies 
of the first-order forces are higher than the ship can 
readily respond to. The second-order forces cause the large 
excursions that must be manually controlled, while the 
first-order forces cause only a ripple on the ship surface 
Srajectory [Ref. 6]. 

The sea also induces moments on the hull of a surface 
Ship. These can be represented by first and second-order 
moments which possess the same basic characteristics as the 
first and second-order forces, and contribute to angular 
Motion in a Similar way. 

The first-order forces exerted by the regular sea have 


the form: 
eae COS (Matto. ) (eqn 4.1) 


wnere 


2. 
We = (Wyaye /9) [(9/yavey > ¥ COsp) | 


The second-order forces FXX, FYY, and moment MZZ have 
constant components and periodic components at twice the 
encounter frequency [Refs. 6, 9, 20]. The second-order 
forces FXX, FYY, and moment MZZ were assumed to approximate 
the constant drifting forces and moments calculated by using 
sea State Program of James L. Cass [Ref. 9] for the Mariner 
Class ship in regular waves. This approximation was: used 
because these were the data available at the time. The 
Sinusoidal part of the the second-order forces and moments 


(frequency of 2 We) was neglected at this point of the work 


17 


[Ref. 6]. The high frequencies of the second-order forces 
and moments are also higher than the ship can respond to. 
They do not effect the ship surface trajectory. 


The exciting forces R and second-order forces for 


i/ 
different encounter frequencies and encounter angle were 
obtained from the sea state program [Ref. 9]. Data input to 
the sea state program for Mariner Class Ship is shown in 
[Rets. al) 10lr 

To see the ship's behaviour in sea state, the simulation 
program was run for ship speed 15 Knots, encounter angle 120 
degree, encounter frequency 0.64 Yradian/Ssecond) 27 -s 
state 8. First and second-order forces and moments were 
added into the surge, sway, and yaw equations that were used 
in the simulation program in Appendix C. Time responses #6r 
u, Ww and § are presented in Figure 4.7, 478) anda sa 
these simulations, optimal gain parameters found by the 
BOXPLX subroutine were used with autopilot. 

The first-order forces and moments are sinusoidal, 
therefore they create only oscillatory ship motions along 
the surge, sway, and yaw axes. The second-order forces and 
moments are constant, thus creating drift motions along 
these axes. 

ie aks apparent that as the encounter frequency 
increases, the effect of the first-order forces on ship 
motions decreases because of the high ship inertia ( acting 
as a low pass filter ) [Ref. 10, 21, 22]. An increasing 
encounter frequency means a decreasing wavelength-to-ship 
length ratio. The constant value of the second-order forces 
increases when the wavelength-to-ship length Laie 
decreases. Regarding the theory, forces and moments are 


significant for short wavelengths.[Refs. 9, 17] 
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Figure 4.7 S.State = 8, Speed=15 Knots, D Type Cont. 
Encounter Angle=120 Deg., Encounter Frequency=0.64 


Yaw Degree vs. Time. 
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Figure 4.9 S.State = 8, Speed=15 Knots, D Type Cont. 


Encounter Angle=120 Deg., Encounter Frequency=0.64 


Speed with Increasing RPM vs. Time. 
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V. COST EUNeI ten 


It is accepted that the performance objective of the 


system in open-sea course keeping is minimum added resis- 


tance due to steering 


For Routine predictions Cilmi. 
fair-weather trial speeds for ships, the increase of resis- 
tance due to steering is often taken to be one percent of 
the towline resistance. Although this may offer little for 
improvement, when it is based on a fuel cost O.5 percent oF 
10 million dollars a year still equals 50000 dou 
Moreover, in moderate to bad weather the benefits of good 
steering rapidly become more evident. It is difficult to 
measure the changes of rate of fuel consumption due to 
different autopilots. A number of alternative cost functions 
have been suggested to describe the rate cf fuel consumption 
[KRefs. 23, 24525). 26562 ecole 

When deriving a cost function for the autoor lou 
requirement is that it must be convenient for ship board 
use. The cost function that has been commonly used in recent 
years is: 


T 
J = lim (472R)ehOw sepa (eqn 5.1) 


‘Ta CO 


> 


This is an approximate form of the exact cost functiem 


Av” EXAGRAEORM 

The performance criterion which characterizes eee 
propulsion losses due to steering may be shown to be that 
derived from excess power consumption per unit distance 
caused by the added resistance due to steering. The added 
resistance due to steering can be related to the surge equa- 
tion where the total instantaneous surge relevant to 


steering is [Ref. 29]: 
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AX = (1/2)X, v2 + (1/2)Xs964 1 (eee ml) Vi 


(1/2 mk ic + Xv (ean 5.2) 


This may be used to determine the energy losses related 
to steering in both calm water and waves. It 1s not conven- 
Zweeror ship board wse, but may be used when finding fuel 
cost and comparing two controllers in fuel consumption. 


Then the exact form of the performance criteria ls: 
T 
ey Z - 2 
J = lim 1790) | 72% saul lal/2 135 pod sO rea) iis 


Nl 


T(E 2a ome en) oo X,svdldt (eqn 5.3) 


The lower limit 1s written as zero, but while using a 
Simulation program, the cost function will be calculated 


after to the transient response time of the system. 


B. APPROXIMATE FORM 

Accurate knowledge of the nonlinear coefficients in the 
exact form Sel —Pocrrotnanee clateriom 1S required for 
Pecwrate results. In addition the criterion itself suffers 
from the disadvantage that the sway velocity is not measur- 
able in practice. 

How can the approximate form which is suitable for ship 
Meard be found {[Refs. 10, 23, 24]? 

Since the sway velocity of the ship is small, the term 
including the square of the sway velocity can be neglected. 
It is also seen that X5v0 and (1/2) (X,,+mX_) r¢ terms are 


small compared to others [Ref. 6]. After these assumptions, 


ZS 


the cost function for the Mariner Class Ship Can Dewy? eae 


as: 
jt 


J = San crazy f Leaargg6" a (A, nm) Vacrmene (eqn 5.4) 
TT © 
@) 
v and r are difficult to measure on ship board. vr can 


be defined in more convenient form for ship board use. 
Yawing and swaying of the vessel is assumed to stem from 
either self oscillations due to the steering system or the 
forced oscillations due to waves. The ship motions resulting 
FEOm oselil later, seaway disturbance are — oscillavueue 
Although these oscillations do not exactly fo0llow a Saas 
soildal pattern, they may be approximated as regular yawing 
of simple peritediiewiee ai Serum rate, sway velocity, and 


rudder angle can be represented as: 


eS See ewan) 


rc 


v Se SBA 3) CTO. ) (eqn Seaa 


Or 
tl 


OF Sin(wttds) 
Theretore: 


J = (1/2) K06* Sin*(wt+ds) 
arin) (Vt ey) COS (Oe ae) (eqn 5 aH 


It is assumed that during low frequency oscillations 
COS Gre Or) aay, i.e. , yaw rate and drift angle are in phase 
With One. 2nG@tneraw ma writing B=v/U as drift rate Wawa 
assuming small amplitude oscillations around the pivot point 


Dp thenjest ron coures sa. 
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(OP/R) = (OP/L)*(L/R) (eqn 5.7) 


Lee) 
iI 


where 


™ 
aT 


(U/r) 


Therefore: 


eaten (OP/L)(L/R) = r° OP (eqn 5.8) 


where 


SS = ity 


Then the expression for the cost function becomes: 


Ak 
ein (1/27) { (01/2) %536? 


Tw 
0 
+(X ptm) OP w*y2 iat (eqn 5.9) 
AWetapor: 
Ay 
J = lim zany f ay?82rat (eqn 5.10) 
T—. © 
6 
where 
(Xo +m) OP yg)” 
se w 
6 
(1/2 )Xgs 
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The accuracy of the ship's path is irrelevant on the 
@eemescea. The important thing 21S minimizing the propulsicn 
power losses. These propulsion power losses can be caused by 
added drag due to steering and external forces and moment 
Piemco Sea waves and wind acting on the Saw. 

Because we and § values are easily measured, ine 


MaenUM COneuo maramecers Canwee Deund by using tne aprrox- 
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[Myace =Crm Of the performance criterion. Howeve 
Mee oreovide evaluation of <xcuel ccnsumpt 
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MpeeeccyYy Speea, and is fwee rrom™iateral disturbance while 
M@magertaking part of a voyage. Moiese aseere ae Ship Ssheuwrd 
Cake an elapsed time t, to cover a distance s, and speed u. 
Meeeccality, the ship does encounter lateral disturbances 
Wnich causes it to yaw at angular velocity r, sway sideways 
With velocity v, and rudder resistance is increased while 
keeping the course. This causes the ship to travel a longer 
Beem (AS) than necessary,a longer time (At) is required 
than for the undisturbed case, and there is speed loss 
e200) - 

The engine useS a certain amount of fuel in a given 
time: 


= W : 
P ie (eqn 6.1) 


E> 
W) 


Fuel Consumption = P G Ceqn 6.2) 
In the disturbed case: 
Fuel Consumption = P <a (eqn 6.3) 


where 


Se 


The ship path can be simulated with the computer simula- 
tion program in Appendix D. And the time to travel the same 
distance along the x-axes can be determined for calm water 
and for the sea state case. 

In fact, to find how much fuel is consumed, the enganw 
specifications must be known. But, when different control- 
lers are compared to determine the percent fuel savings 


in the same situations, there 1s no need for engine 


specifications. 
> Hue Peon ctiijet ton Pt") eeeees 
Rawr SS a See x 100 (ean uae 
{ 
Pt 1 
where 
Pet', = Fuel consumption at t',, elapsed time "toute aes 


the given distance along the x-axes while using controller 


#1(assume as an automatic optimal controller) 


Pet's, = Fuel consumption at t 5, elapsed time tour eawem 
the given distance along x-axes while uSing controller 


#Z2(assume as a PID controller) 
them: 


% Euel Consumption cle a ule 
eo SS ee x 100 (ean 6.5) 


si) 


The sea state forces and moments, and added drag due to 
steering affect the’ elapsed time to travel the given 
distance along the x-axes. Mice ratiQu.giyves the tfuel 
consumption ratio including all sea state forces, moments, 
and added drag due to steering. But, the effect of added 
drag due to steering for different controllers is needed to 
compare fuel consumption. 

The other disadvantage is that it can be used only for 
small sea state disturbances, since the ship speed decreases 
too much in high sea state. It may be even sufficient to 
drive the ship backwards [Ref. 6]. The wilecil (Dice Seis lene Ticlel= 
decreasing speed is that this changes the values of the 
hydrodynamic coefficients of the nonlinear equations given 
in Table 2, 3, and 4 in Appendix B. The big speed differ- 
ences from steady-state speed, 15 Knots, makes the nonlinear 
coefficient values incorrect to represent the ship-steering 
dynamics. 

Pome ees me tor small disturbances, i.e. , when the 
representation of ship-steering dynamics does not change, 
then this method may be used. But, for large sea distur- 
bances, the ship speed must not be allowed to decrease much 
to affect the representation of ship-steering dynamics in 


the simulation program. 


B. CONSTANT SPEED 
Using Table 6.1 [Ref. 6], it is seen that the added drag 
force due to increasing RPM is a linear function. The linear 


" added force vs. RPM " is found, and used in the 


equation 
Simulation program to Keep the speed constant. 

While attempting to maintain speed constant, it is more 
werercult to find the fuel consumption ratio, since P is 
Memying. P cannot be factored out, as in Eqn. 6.5. 

How can a method to find the fuel consumption ratio 


Without using engine specifications be found? 


3k 


The added resistance due to steering is: 
= Z 2 7 
AX = 7 + (A 2 ee t (1/2 )Xes6 


+ (Xyptm)vr + X evo (eqn 6.6) 


The total surge equation is: 


“total ~ “ * “ealm (eqn 6.79 


where 


Xx =T (1-=t) 


Gaim 


The fuel consumption ratio of added resistance duemigs 


steering to the total surge equation, assuming Comoe 


overall propulsive efficiency, [Refs. 2, 25] is: 
% Fuel Consumption Ax 
pei - a x 100 (eqn 6.8) 
soul 


Using this ratio, it is possible to compare the fuel 
consumption of the Automatic Optimal Controller and jwise 
that of the PID. To find exact fuel consumption, engine 
specifications are still needed. This ratio will be used in 
Chapter 7 to compare the PID with the Automatic Optimal 


Conmtro) ber. 


a2 


Shes Ib) 3, al 
SPEED AND eDED SORCE VS. RPM 


RPM Speed Added Force 
(Knots) (LB) 
SS a 133140 
S/S) ts | 9 | SS 70 
Beet | oe | Z21Ome 
294s | LOz.O | ZO 
Sila | tO) | 31400 
628 zeae. | os 7380 
64.6 | ipsa @ | 45500 
Sons | 14.0 | 335) (S10) 
68.6 digs ol) 63420 
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VII. COMPARISON OF PID AND AUTOMATIC OPMIMAL CONTROL s 


The Automatic Optimal Controller offers the potential 
for minimization of propulsion losses due to steering in the 
open sea, and for removal of operator judgement from 
steering control of the ship. The optimal parameters which 
were found by using the BOXPLX subroutine in the Op@iiiam 
Gain Program were shown in Table II [Ref. 1]. The approxi- 
mate form of the cost function was used in the ~BOQxaims 
subroutine. 

Time responses of wy and § are shown in £igureuam 
through figure 7.10 for different encounter angles using the 
optimal parameter values in Table II 

No provision for automatic adaptivity to either specs 
load or seaway exists in the Universal Gyropilot, PID. Some 
adjustments of control parameters are possible, however, 
through an operator interface. Details of such adjustmenue 
and the structure and parameter values of a PID controller 
are contained in Appendix A. 

Although a PID controller can come close to emulating 
this performance under specified internal and operator 
settings, the operator has no means’ to judge steering 
performance other than by observation of the course 
recorder. The normal tendency of the operator 1s to change 
the external controls available to him to bring about the 
reduction ian heading “erreue It greatly increases the 
propulsion losses due to steering. 

The goal in this chapter is to compare an Automatic 
Optimal Controller with a well tuned PID controller in fuel 
Savings.Now the question is how the Well-Tuned PID can be 
defined. There are two ways to answer this question. One way 
is to define a well tuned PID as an Optimal PID Controller 
and to use the same performance criterion used in finding 


optimal parameters of a D Type Automatic Controller, i.e., 
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the minimization of added resistance due to steering. The 
PeeimeaecOnrLOolL Pabanecers OL a PID in different sea state, 
speed, encounter angle, and encounter frequencies were found 
by Using the simulation program to compute the cost function 
of all combinations of internal and external control parame- 
ters, making the modifications in the simulation program 
shown in Appendix C. The results were shown in Table III 

em seo nOnsoomOrmuemancdeg are shown .in Figure 7.11 
through 7.20 for different encounter angles using the 
optimal parameter values in Table III . As it is seen from 
the figures, the Automatic Optimal Controller and PID with 
optimal parameters after the transient time makes the ship 
yaw sSinusoidally around the commanded course. To see their 
effects in added resistance due to steering, and in fuel 
Savings, the comparison of them was shown in Table IV . The 
added resistance due to steering and fuel consumption ratio 
are almost the same, since the same performance criterion 
was used, and all internal and external controls of PID are 
adjusted, not only the external controls. 

The PID with optimum parameters cannot be used as a well 
tuned PID, since the operator could not find the optimum 
parameters by himself to minimize the added resistance due 
Pemsteering and also could not change the internal control 
Ieaiamecers. But it is clear that if the optimum control 
parameters are found in the same way as for the D Type 
Controller, it gives the same fuel Savings as an optimal D 
Type controller. 

The other way is to approach the question as an oper- 
ator. There are no means to adjust the external controls 
available to the operator to minimize the propulsion losses 
due to steering other than by observation of the course 
recorder. The operator can adjust the external controls to 
reduce the heading error, the only system areegenne available 
to him. The external controls available to the operator are 


Weather Adjust Gain, Rudder Multiplier Gain, and Rudder Rate 
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Multiplier. To reduce the heading error, the Weather Adjieg 
Gain is set to "0", and to get the fast response of rudder, 
the Rudder Multiplier Gain is at "3", and the Rudder rate 
Multiplier setting is “1”. These correspond (To Ki 1s) a 
is "3" “and the =@valitemern T, is half of the internal comemem 


setting. 


The 21nternalmieome sous Th, T2, and T3 are not known for 
the Mariner Class Type ship. These internal control parame- 
ters are found by running the simulation program used in 
finding optimum parameters for the PID with K, and Ky set to 
"1", and the proportional part of the PID is omitted in the 
control law of PID, Egn. A.1. The results are following. 


n= 25 
My SL 
hy = eloyo 


Time responses of yw and § are shown in Figures 
Cheeugie 7 aoeF As 1s seen in Figure 7.22, /.24, 7.26, 7ieaee 
and 7.30, the well tuned PID uses a larger rudder angle to 
reduce the heading error quickly. The rudder angle alse 
oscillates in large magnitudes. It tries to get the ship to 
the command heading angle as fast as possible. It increases 
Ene WeCcee So consequently the fuel cost and added resis- 
tance due to steering increase. The comparison in fuel 
Savings and in added resistance due to steering 1s shown in 
Table V . Although the optimal internal control paranere = 
are used in a well tuned PID, there are differences in fuel 
consumption ratio and in added resistance due to steering, 
especially at quartering sea waves, i.e., 030-060 and 
120-150 degrees. The Automatic Optimal Controller provides 
fuel savings in excess of 0.5 % over a well tuned PID when 
operating at the design speed at random headings in sea 


states. 


S16 


Hiemeeieehna, COmterol parameters of a well tuned PID in 
the deterministic model are represented by optimal values. 
They might not be the actual values used for a Mariner Class 
Ship. The differences in fuel consumption ratio and in added 
resistance due to steering should increase if actual values 
are not equal to the optimal internal control parameter 
values used in the simulation computer program. It is also 
possible that the different operators can set the Rudder 
Multiplier Gain and Rudder Rate Multiplier in different 
ways. The fuel consumption ratio and added resistance due to 


steering also changes in relation to these valtes. 


Sy 
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TABLE IV 


COMPARISON OF OPTIMUM PID WITH AUTOMATIC 


CONTROLLER 


Sea State = 8 

Ship Seeec = 15 Knots 
Waviemerequency = 0.53 
Added Resistance (libre) 


Encounter “Angle 030 deg. 


D Type Spcimum PID 
Added Resistance 
Due to Steering 5547 5892 
Fuel Consumption 
Ratio ( % ) Sie IONS. Ses y4 


Encounter Angle O60 deg. 


D Type Optimum PID 
Added Resistance 
Due to Steering 52422 92734 
Fuel Consumption 
Ratio ( Son oa] 34.080 


Encounter Angle O90 deg. 


D Type Optimum PID 
Added Resistance 
Due to Steering 108930 109250 
Fuel Consumption 
Ratio ( % ) 59.749 59.928 


Encounter Angle 120 deg. 


D Type Optimum PID 
Added Resistance 
Due to Steering ps7 7 Sulesiees 
Fuel Consumption 
Reawemo ( + ae, 179 28.442 


Encounter Angle 150 deg. 


D Type Optimum PID 
Added Resistance 
Due to Steering 14756 ie 21 
Fuel Consumption 
Ratio ( ¥ ) Se OsZ 2) 5 SHON) 


CRE, 


OPTIMAL 


Difference 


She egeys Vi 


COMPARISON OF WELL TUNED PID WITH AUTOMATIC OPTIMAL 
CONTROLLER 


Sea State = 8 
Siovalje Sh = 15 Knots 
Wave Frequency = 0.53 


Encounter Angle 


Dp elyipe 
Cént 
Added Resistance 
Due to Steering 5547 
Fuel Consumption 
Ratio ( % 2) 6 NOs 
Encounter Angle 
D Type 
aoa) 
Added Resistance 
Due to Steering 52422 
Fuel Consumption 
Ratio ( 7 33.947 
Encounter Angle 
D Type 
Sone: 
Added Resistance 
Due to Steering rOe230 
Fuel Consumption 
Ratio ( & ) 59.749 
Encounter Angle 
D Type 
ee 
Added Resistance 
Due to Steering Sls 77 
Fuel Consumption 
Ratio 7 ) Fuos | ft, 
Encounter Angle 
D Type 
Sent 
Added Resistance 
Due to Steering 14756 
Fuel Consumption 
Ratwo (a7 BS. 052 
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O30 deg. 


Well Tuned — 
PID 


O60 deg. 


090 deg. 


59.934 

120 deg. 

Well Tuned 
PID 


Difference 
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Figures 7a S. State=8, D Type Cont., Specd—15e nee. 


Encounter Angle=060, Encounter Frequency=0.43 


Rudder vs. Time. 
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Rudder vs. Time. 
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Figure 7.8 
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Beegure 7.10 Se oeare—o, Di Type Cont., Speed=15 knots 
Encounter Angle=150, Encounter Frequency=0.67 


Rudder vs. Time. 
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Rudder vs. Time. 
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Figure 7.14 S. State=8, Optimal PID, Speed=15 Knots 
Encounter Angle=060, Encounter Frequency=0. 43 


Rudder vs. Time. 
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Figure) 7220 S. State=8, Optimal PID, Speed=15 Knots 
Encounter Angle=150, Encounter Frequency=0.67 


Rudder vs. Time. 
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Figure 7.22 S. State=8, Well Tuned PID, Speed=15 Knots 
Encounter Angle=030, Encounter Frequency=0.33 
Rudder vs. Time. 
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E re 7.24 S. State=8, Well Tuned PID, Speed=15 Knots 
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Figure /s26 S. State=8, Well Tuned PID, Speed=15 Knots 
Encounter Angle=120, Encounter Frequency=0.64 


Rudder vs. Time. 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSTONS 

An A Type controller was used in finding optimal parame- 
ters to minimize the propulsion losses due to steering in 
previous theses [Refs. 8, 10, 21, 22]. But to prevent the 
drifting effect of second-order forces and moments on the 
ship's motions, a D Type controller was used in designing 
the Automatic Optimal Controller. The only difference 
between them is the integrator part omitted in the A Type: 

Savings in fuel for an Automatic Optimal Controller yaa 
Optimal PID, and a Well-Tuned PID were compared by uSing a 
new approach explained in Chapter 6 that does not use engine 
specifications. It can be used in finding any ship's fuel 
Savings by finding the added resistance due to steering. 

The optimal internal and external settings of a PID to 
minimize the propulsion losses due to steering are found by 
the modified simulation program shown in Appendix C, since 
there is not enough information about the settings .of 
control parameters of PID for the Mariner Class ship. The 
Same performance criterion explained in Chapter 5 was used 
in finding the optimal parameters for both of them. As a 
result, it is found that the optimal PID iS as: good as 
Automatic Optimal Controller. They provide almost the same 
fuel savings in the same situations. 

The well tuned PID was defined as a PID controller to 
reduce the heading error as fast as possible and to get fast 
response of the rudder to do this. Only external control 
settings (explained in Appendix A) are available to the 
operator. The optimal internal control settings were used in 
Simulations.The Automatic Optimal Controller provides fuel 
Savings in excess of 0.5 % over a well tuned PID when oper- 
ating at sea state, 8, speed, 15 Knots, and wave frequency, 


O.53 radian/seconds. 
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Be RECOMMENDATIONS FOR FUTURE STUDIES 

Additional work should be done to compare the Adaptive 
Automatic Optimal Controller with a well tuned PID for 
different sea states, speeds, and encounter frequencies. 

In this thesis the well tuned PID was used as explained 
in Chapter 7. If more information is available for the PID, 
especially if a precise definition can be found for the well 
tuned PID, then the comparison in fuel savings with these 
well tuned PID settings will provide better insight into the 
value of the Optimal Controller. 

Li ties unesis, SUage, Sway, and @yaw equations of motion 
were used in simulations. The roll equation of motion may 
Bemeccead to the simulation program and the optimal parame-— 
ters of the Automatic Optimal Controller may be found by 
using the BOXPLX minimization subroutine. Then the compar- 


ison in fuel savings will be more realistic. 


2) 


APPENDIX A 
PID CONTROLLER 


The standard form of a PID controller is baseauie 
proportional-plus-integral-plus-derivative contiou= No 
provision for automatic adaptivity of the control lorie 
speed or seaway effects exists in the PID. It relies on 


adjustment of the control parameters by the operator. 
1. GENERAL DESCRIPTION - 


An autopilot commonly found aboard merchant ships has 
the characteristics of a PID controller. The PID control 
is described by [Refs. 23, 24]: 


Glsies)\- I'3s 


where the gains are determined by operator interaction in 
the following manner: 
1. Weather Adjust Gain 
K, = 1/3 within "weather wadyus = wom 
= 1 outside "weather adjust" zone 
The weather adjust zone is variable from 0.3 to 5.0 
degree heading error. Panel is markec ©) seems = 
2. Rudder Muleiplier Gama 
Ko = 1 to 3, continuously adjustable (panel 


COnELomD 


SNe. 


Ss intemnaleconkrol Settings 
The time constants Oey V5, and T3 are preset for the 


specific ship, with the exception that may be reduced by a 
factor of 2 by operator interaction. The ranges of possible 


gyropilot time constants are: 


1. Derivative Time Constant 
1) See Om z ope, OO 150, 200 second 
"Rate Mult" control reduces the time constant by a 


factor of 2 when at minimum setting. 


2. Derivative Filter Time Constant 
Ts Soa eo S,. LOpe lS second 


3. Integral Time Constant 
Tz = 200, 400, 600, 1000 second 


fee APPLIED TO SIMULATION PROGRAM 


Figure A.1l corresponds to controller PID which has the 
morm {[ref. 22]: 


(1+T5s)4 Tas 


This can be written in simulation program as this: 


DX2 = ( YAWE - X2 ) / Tp 
ORG = (OD ea? oe 


7“ 


Where 


<2 ee) + D2 ELT *} 
X3 = X38 + ( DX3 * DELT ) 
X4 = X4 + ( YAWE * (K,/T3) ) * DELT 


SNe, 





PID Controller in Simulation se regea 


Figure A.1l 
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APPENDIX B 
NONLINEAR EQUATIONS OF MOTION AND COEFFICIENTS TABLES 


Nomemeiear equacions of Motions: 


eect > Xe yu = Ea (u,Vv,r,o) 
manag. : (m = Ys)V fp eee ee A 8G: Cue eo) (eqn B.1) 


‘ 
Z, 
ry 

H 


N-Eq.: (mXq -Ne)v + (1 Bere; vO) 


Where 


S 


mv, 2,0)— X* + X¥, Au + (1/2)X,,,Au* + (1/6)X a> + 


uuud 
: 2 2 2 

BO ye | C72) Xpemxe ye PMGI/2)X5n5* + 

U2) en + (1/2) eee AU (1/2)Xa06“Au + 

Con) ve X., sv + X5rd + X 


vruvFAu + 


Xs, voAu tr Xs yrdAu 


ey 


fo(u,v,r,0)= Y* + ¥*,Au + ¥*,,Au> + Yow oe 
(1/2 eee (1/2) Yygsvd" + Yo vAu + 


(1/2)¥y ,vAu2 + (Y,-mu)s + (176) 


(1/2) eee 1 (1/2)¥,99rd* 2 Y ey na + 


(1,2) ¥puyrdu® + Yd + (17/6) ¥ 535° + 
(2) oe Se omel (1/2 ea one 


Ya aol = OPE 


v3 
vvv" 


£3 (U,V 70 eee ee eee N* Au? + No lean + 


(1/2)Njppvr* + (1/2)Nygevd* + Ny vAu + 


(1/2)Nu yy vAue # (N-mX ue en eee 
(1/2)Npyyrv? + (1/2)Niser§* + NirAu + 
(1/2)N,.,rAu* + Nsd + (1/6 )Ngasd? ‘ 


(1/2)No.6v% + Ng SAu + (1/2)Ng SAU + 


N-gvrd + (1/2 )NopOr4 


e172 


ay Lor 
Expansion. 
and Dynamic 
response 
Terms 


(m-X-) 


xy 
(DY 


(1/6) Xiy 


(ae 
(1/2 5 sue 


(1/6) Xge 


ye eae 
(1/2) Xp yy 
(1/2)Xee, 


(X..+m) 
Xs 
Xs 


Ky, ru 


AS u 


rou 


xX* 


TABLE VI 


COEFFICIENTS OF X-EQUATION 


Paenei fier 


Blk gl 


Fortran 
Program 


cil 


en 
X2 
X3 


X4 


oS 
X6 


X7 
X8 


XO 
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Nondim. 
Paccor 


(1/2)p 


(1/2)p 
(1/2)p 
(1/2)p 


(1/2 )p 
(172) p 
(1/2)p 


(1/2)p 
(1/2)p 
(1/2)p 


(1/2)p 
(1/2)p 
ey 20 
(1/2)p 
(1/2)p 
(1/2)p 


(1/2)p 


LBP? 


LBP“S- 
LBP~ 
LBP“ /S 


LBP# 
LBP* 
LBP“s¢ 


LBP“ /S 
LBP*/s 
LBP“S 


LBP? 
LBP“s 
LBP?S 
LBP» /S 
LBP“ 
LBP? 


LBP“s¢ 


Nondim. 
Goctt. 


*10° 


840.0 


= JAC aI8) 
45.0 
m© 3 


= eSiSi cas. 
1S 6. 
-94.8 


ESS @ 
So 


Taylor 


Expansion. 
and Dynamic 


response 
Terms 


Cie) 
(mans) 


“SG 


76 
(1/2) Y ory 
(1/2) a8 


(Y,-mu ) 

WIE) Severs 
(172) ee 
(1/2)Y,5¢ 


*§ 
(1/6) Y 550 
(1/2)Ys5.4, 
(1/2)Y¥srr 
¥8u 

(1/2) Ysuu 


Yyrg 
yx 


Nee 


Jefe ley gue al je th(ey 1s 


TABLE VII 
COEFFICIENTS OF bourse) 


Bh @) 


Fortran 
Program 


es 


ou 
eZ 


x3 


Y4 


bas: 


YO 


SCG 
me 


xO 
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Nondim. 
Factcon 


(1/2)p 
(1/2 )p 


(1/2)p 
(1/2)p 
(1/2)p 
(ea 


(1/2)p 
(1/2)p 
(1/2)p 
(1/2)p 


(1/2)p 
(1/2)p 
(17 Zp 
(1/2)p 
(1/2)p 
(1/2)p 


(1/2)p 


(1/2)p 
(1/2)p 


LBP* 
LBP@ 


LBP“S 
LBP“ /s 
LBP*/s 
LBP“sS 


LBP~s 
LBP~/S 
LBP? /S 
LBP°sS 


LBP“s? 
LBP“ S$ 
LBP# 
LBP* 
LBP“s 
LBP# 


LBP? 


LBP“s¢ 
LBP“Ss 


Nondime 
Coeff. 


*10> 


1546. 


“Oe 


-L1lG6Gr 
-807 ce 


aS) 


-499. 


153567 


Zige 


2 2)0). 
S'S) - 


38 
OF 


(On 2 SS. 


QO OV OO 


MO O ND -& 


6 
0 


PAS bee 1 1 
COmmnre PENtTS. OF N=EQUATION 


Taylor. Tdentifier Nondim. 
Bapals LOM ay Ba Gic@ Ie 

and Dynamic Fortran 

response Program 

Terms 

(mX,-Ne ) C5 (1/2)p LBP* 
(I, -Ne.) C4 (1/2)p LBP? 
N,, N1 (W772) paEBEs 
(1/6)N N2 (1/2)p LBP?/s 
AyZ Ne (1/2)p LBP?/S 
(1/2) Ny gs N3 (1/2)p LBP?S 
(N,.-mXgu) N4 (1/2)p LBP*s 
(1/6)Ny ny (1/2)p LBP®/s 
@y72N N5 (1/2)p LBP*/s 
(1/2 )Nyas (1/2)p LBP*S 
Ne N6 (1/2)p LBP?s4 
(1/6 ) Nags N7 Gy 2 penBe s- 
(1/2) Ng vn, N8 (1/2)p LBP? 
(1/2 Nery (1/2)p LBPts 
Ngu (1/2)p LBP?s 
(1/2)N vs (1/2)p LBP* 
N* NO Ul/2ya Wiis | 
Ne, (1/2)p LBPYs 

| 
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Nondim. 
Goer. 


*10> 


SAV aed) 


SZ. 


S| oye 
neo 6. 


Le 


OG. 


-5483. 


gad oo 


Sa on 


Cy @ ie ® By &. 0) i © ee en od 


Oy (SO AS I 28 9, 


Ke) 
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APPENDIX C 
SIMULATION PROGRAM 


Cx xx DECLERATIONS OF VARIABLES KRRREKRKKKKRRKRRKKKRE KR SE eee 
REAL MH(10000),1T(10000) 
REAL*8 L,L2,L3,L4,L5,L6 —— 
REAL*8 X,XDOT,Y,YDOT,U,UDOT,V, VDOT, YAW, R, RDOT 
REAL*8 X1,X2,X3,X4,X5,X6,X7,X8 
REAL*8 YO,Y1,Y2,Y3,Y4,Y5,Y6,Y7,Y8 
REAL*8 NO,N1,N2,N3,N4,N5,N6,N7,N8, TEM, TEMP, TEM1 


C TEMP 1 
REAL*8 C1,C2,C3,C4,C5.F1,F2,F3,F2x, F2Y Nez) Ea 
C /MZZ | 


REAL*8 RO,DELT,S,DU,U1,K,Z1, 2252162, KD De 
REAL*8 DYAWE,YAWE, YAWC,1ISR,ISE,TDIFF, LAMDA,RPM,RES 
GS , EHP 

REAL*8 S$1,$2,DS$1,DS2,D, LOST, YAWDEG, DDEG, $3 

REAL*8 MASS,1IZ,XG, YVDOT, NVDOT, YR, YRDOT, NR, NRDOT, FX 


c Ee 
REAL*8 RX,RY,RZ,TX,TY,1I2,WA, WE, RXR, RYRPRA TD RY eee 
S plz 


REAL*8 MACGMAX,MACMIN, TIME, ETIME, DELST1, DELSITZ 3 bens 
REAL*8 BARFUE, LBFUEL 
REAL*8 TOTFRI, FOTERZ, PERCI fone 

C***WHEN PID CONTROLLER IS USED , ADD THE FOLLOWING 


C***xVARIABLES Hi i HK 
c REAL*8 PIDV2,PIDV3,PIDV4, PIDDV2,PIDDV3, PIDT1, PIDT2 
Cc C /pIDT3, PIDK1, pidk2 


C***WHEN OPTIMAL PID PARAMETERS ARE WANTED TG EIND? 
C***ADD THE FOLLOWING VARIABLES AND GIVEN BATA. VAGUE. 


e REAL*8 KK1(2),KK2(5),TL1(9) ct sn eee) 
6 DATA Kis) ©ascy la, 
ec DATA KK2(1),KK2(2),KK2(3),KK2(4) , KK2( 57 Eee 
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é Py cy cme 
Poe een 2) Tr 3), TT1(4),1TT1(5),TT1(6) 

€ PTT ( 7m TT1(8), 

Cc Brio (Sess oe IsaoS.,50.,100.,150.,200./ 
wee Tn (1) T6272), TT2(3) , TT2(4) ,TT2(5)/1.5,2.5,5- 

eS moO. , 152 
Doub wns) Tra(2),5bo (3), 1T1T3(4)/200. , 400. , 600. 

é | OO. / 


CEKKRKRKREKREKRKEKKRRRKREKRRREKKKRRRKRERRKEKRRKRRERERKERERRERKEREREREKERKRREE 


© ae @ A es nes: PD pe Ga Se 


c 

C . 

Se Tahoe Plooer LA STATEMENTS ARE USED ONLY WHEN *****% 
oS RUNNEING TO CET GRAPHS. THEY ARE NOT USED oe * * 
C*** WHEN FINDING OPTIMAL PARAMETERS FOR PID. BES eg 


CALL TEK618 

Camm PAGE(11.0,8.5) 

CALL NOBRDR 

CALL HEIGHT(0.15) 

CALL MX1ALF('STANDARD','&' ) 

CAM MX2ALF('L/CSTD','#' ) 
CERKRKRKEKRKRRKERRRREKRRKERRERRKERERKEKREKKRKEKRRRRKRREKKEKKRKREKEEE 
S 
C*** USE THIS PART TO FIND OPTIMAL PID PARAMETERS ****** 
DO 11 IK1=1,2 | 
Domes IkK2—1 5 
DO 3S" 1T1=a, 9 
Demet r2—1,5 
Domes IT3=1,4 


(cal @ ND ais @ aes ia @ a 


C*** WHEN FINDING OPTIMAL PID PARAMETERS, *****# 
C*** USED ETIME=1000. Re RE 
ETIME=2000. | 
TIME=O. 
ITIME=0 
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ICOUNT=1 


C*** THESE ARE OPTIMAL PARAMETERS OF AUTOMATIC@G3@P Ti) 17 [aa 
C*** CONTROLLER FOUND BY RUNNING OPTIMAL GAIN PROGRAM a 
C*** IN REFERANS 7, 27021522 ee ce 
K=0.6983 
Z1=0.0102 
PI=34- 3620 
KD=569.1469 
C*** HERE, YOU CAN USE SEER 20) eee PARAMETERS OR “5 


C*** OPTIMAL PID PARAMETERS FOUND BY Sau ee ia Bie os SS 
C*** SIMULATION PROGRAM USING THE PAR eons ene |e 1 
C*#** OPTIMAL Pise yArtUESe bless 
Ss PIDK1=1.0 

C PIDK2=s30 

C PIDT1=25. 

C PiDTZ=l>- 

Cc PIDT3=200. 

S 


C*** INITIALIZE ALL THE NECESSARY VARIABLES ******% eee ke KH 
C INITIALIZE THE GOstr HUNemnen 
ISE=0.0 
ISR=0.0 
sGOLIAIT =) © 
LAMDA=5 . 48 
C X,XDOT,Y,YDOT ARE FIX COORDINATES ON EARTH 
KONO 
V=Om0 
ED CHE 0), © 
VDOT=O0n0 
C U,UDOT,V,VDOT ARE FIX COORDINATES ON SHIP 
WO 1c 
UDOT=0150 
VDOT=0.0 
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Q © 


On 
F2=0.0 
F3=0.0 
YAW=0.0 
YAWE=0.0 
R=0.0 
RDOT=0 .0 
RPM=68. 6 
PEo=O. 
EHP=0.0 
BARFUE = 0.0 
LBFUEL = 0.0 
RO=1.9876*.5 
Ca92 1! 
YAWE=0 .0 
TEM=0.0 
TEMP=0.0 
TEM1=0.0 
TEMP 1=0.0 
Dol=Oo 
DS2=0.0 
S1= 0) (6 
S2=0.0 
PIDDV2=0.0 
PIDDV3=0.0 
PIDV2=0.0 
PLD —O20 
PIDV4=0.0 
SMALL=99999999.9 
D = RUDDER ANGLE 
D=0.0/57.296 
ORDERED SPEED IN FEET/SEC 
15.*1.689 FT/SEC=15 KNOTS 
U1=15.*1.689 
AT STEADY STATE ACTUAL SPEED (U} 
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COMMAND SPEED (UC) 


U=U1 

Eb sZ29- 
L2=h=2 
Ne Ye Ur) E> 
L4=L*L3 
L5=L*L4 
L6=L*L5 


Cc FORCES IN X,Y DIRECTION COMPUTED IN eer eas 
C MOMENTS IN 24 

FX=0. 

FY=0. 

M2Z=OQ. 
& 
© . 
C*** SEA STATE FORCES AND MOMENTS FOUND SY RUNG. 
THE JIM CASS'S SEA STATE PROGRAM IN REFERANS 8. SN 
C*** “BIRST ORDER FORCES Nb eae 


Cex 


S 


RXR=-.12887D4 

RAL==. 9051205 

RYR= .74476D3 

Ryi= .1T226bs6 

MAR =e OC 

MZ0—-7, 1OOZED7 

SECOND ORDER FORCES "AND Wi Ciata 
k= SO) SS BSS 

P2y7=". LoeZzen 4 

M2Z= .96895D-4 


C*** FINDING FORCES AND MOMENT MAGNITUDES AND PHASES) a5 


RX=(RXR**2+RXI**2)**_5 
RY=(RYR**2+RYI**2)**_5 
RZ=(MZR**2+MZ1I**2)** 5 
TX=DATAN2 (RXI,RXR) 
TY=DATAN2 (RYI, RYR) 
TZ=DATAN2 (MZ1I,MZR) 
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C 
C 


PXM=Z. *RO*GC*L*AWA**2*EF2X 
Reva RONG a WwA**2*h 2 Y 
Mga 2 ee noe («42 )~(WA**2 )*M2Z 


C*** SEA STATE AND SHIP SPECIFICATIONS ******** ek RRR RKKK KX 


C 
C 


Os @ aes oan aed ae 


Se Pet eColteeWavgedslGHl:-osteeslTATE G-19 FEET, 7-15 
Slee / 

WA=19. 

ENCOUNTER FREQUENCY: (WHEN ENCOUNTER ANGLE IS 30) 

WE=0 . 33 

ADDED MASS AND ADDED INERTIA TERMS: 

MASS=.11155D+07 

bZ2—0 7 0@08 92 ~RO- LS 

XG=-11.96 


200 CONTINUE 


INPUT YAW COMMAND 
YAWC= 0.0/57.296 


TO PREVENT THE SHIP TURNING MORE THAN ONE DEGREE 
PER SECOND AT THE BEGINNING OF THE PROGRAM SINCE 
THE SHIP CANNOT BE IN HIGHER SEA STATE INSTANTLY. 
THERE MUST BE TIME PERIOD C FROM GOING SEA STATE 
ONE TO SEA STATE EIGHT. 
TEM1=YAW-TEMP1 
IF(TEM1.GT.1.0/57.296) YAW=TEMP1+1.0/57.296 
IF(TEM1.LT.-1.0/57.296) YAW=TEMP1-1.0/57.296 
TEMP 1=YAW 


ERROR SIGNAL TO DRIVE RUDDER( YAW ACTUAL - YAW ORDERED) 
YAWE=YAW - YAWC 


S=((U*U)+(V*¥V) )**.5 


gal 


DU=U-U1 
C***IF D TYPE AUTOMATIC OPTIMAL CONTROLLER IS USED, ***** 
C***USE THESE ee ae 
DS1=(YAWE - $1)/P1 
D=(S1 + DS1*Z1)*K+S2 


e 

C***IF OPTIMAL OR WELL TUNED PID IS USED TO GET GRAPHS*** 

Cx* USE THIS PART. Ka 

C PIDDV2=(YAWE - PIDV2)/ PIDT2 

C PIDDV3= (PIDV2- PIDV3 )} 7 PIDTZ 

é D = (PIDK1*PIDK2*PIDT1*PIDDV3) + (PIDK1*PIDK2*YAWE) 

C é + PIDV4 

@ 

C***WHEN FINDING OPTIMAL PID PARAMETERS, USE THIS PART.**** 

C PIDDV2=(YAWE - PIDV2)/ TT2(1T2) 

@ PIDDV3= (PIDV2- PIDV3 ) / TT2(1T2; 

é D = (KK1(1)*KK2(IK2)*TT1(IT1)*PIDDV3) + PIDV4 

é C +(KK1(1)*KK2(IK2)*YAWE) 

é 

C*** THE RUDDER ANGLE COULD NOT BE CHANGED MORE THAN **** 

C*** 2.5 DEGREES PER SECOND. Aw KK 
TEM=D- TEMP 


IF(TEM.GT.2.5/57.296) D=TEMP+2.5/57.296 
IF(TEM.LT.-2.5/57.296) D=TEMP-2.5/57.296 


IF ( D.GT.(35./57.296) ) | D =" 65m oy eee 
IF ( D.LT.(-35./57.296) ) D =eee5eyo7eece 
TEMP=D 
C 
C*** THE VALUES OF HYDRODYNAMIC COEFFICIENTS TO REPRESENT** 
C*** SHIP STEERING DYNAMICS. ( REF. 1, 6 ) i 


C AXIAL FORCE HYDRODYNAMIC COEFFICIENTS (SURGE) 
X1=(-0.00120)*(RO*L2*S) 
X2=(0.00045)*(RO*L2) 
X3=(-0.000103)*(RO*L2/S) 

X4=(-0.00898) *(RO*L2) 


IZ 


X5=(0.00018)*(RO*L4) 
X6=(-0.000948) *(RO*L2*S*S) 
X7=(0.00798)*(RO*L3) 
X8=(0.000832)*(RO*L2*S) 

C LATERAL FORCE HYDRODYNAMIC COEFFICIENTS (SWAY) 

C YO=(-0.000036) *(RO*L2*S*S) 
YO=0.0 
Y1=(-0.011604) *(RO*L2*S) 
Y2=( -0.08078) *(RO*L2/S) 
Y3=(-0.000038) *(RO*L2*S) 
Y4=(-0.00499)*(RO*L3*S) 
Y5=(0.15356)*(RO*L3/S) 
Y6=(0.002779) *(RO*L2*S*S) 
Y7=(-0.0009)*(RO*L2*S*S) 
Y8=(0.011896)*(RO*L2) 

C MOMENT ABOUT Z-AXIS HYDRODYNAMIC COEFFICIENTS (YAW) 

C NO=(0.000028)*(RO*L3*S*S) 
NO=0.0 
N1=(-0.002635)*(RO*L3*S) 
N2=(0.016361)*(RO*L3/S) 
N3=(0.000125)*(RO*L3*S) 
N4=(-0.00166)*(RO*L4*S) 
N5=(-0.05483)*(RO*L4/S) 

Miic—(-OsO0dces + (RO«ES*S<S) 
N7=(0.00045)*(RO*L3*S*S) 
N8=(-0.00489)*(RO*L3) 

C COMMON COEFFICIENTS 
C1=(0.0084)* (RO*L3) 
C2=(0.01546) *(RO*L3) 
C3=(-0.000086)* (RO*L4) 
C4=(0.000829)*(RO*LS) 
C5=(-0.000227)*(RO*L4) 


Cc REGULAR WAVES 


73 


OQ 


Ck*e* 
Ck*k* 
Ck*eEK 


Ck*x* 


GC 
C 
Cxx* 


Ofojoit 


EX=WA*RX*DCOS (WE*TIME+TX) 
FY=WA* RY*DCOS(WE*TIME+TY) 
MZ=WA*RZ*DCOS (WE*TIME+TZ) 


IF(TIME.EQ.0.0)FX=0.0 
IF(DABS(FY).LT.0.00000001) FY=0.0 
IF (DABS(MZ) .LT.0.00000001)MZ=0.0 


EQUATIONS OF VMGETe! 


Fl = X1*DU + X2*DU*DU + X3*DU*DU~DU + X4-*\ a eee 


1 + X6*D*D + X7*V*R + X8*V*D Role ees 
F2 = YO + Y1*V + Y2*V*V*V + Y3*V*D*D + Y4*R 
+ YS*R*V*V + Y6*D + Y7*D*D*D + Y8*D*V*V 
+ FY + FYY | 
F3 = NO + N1*V + N2*V*V*V + N3*V*D*D + N4*R 
+ NS*R*V8V + N6*D + N7*D*D*D + N8*D*V¥V 
+ MZ + MZZ 


ADDED RESISTANCE RELEVANT DUE STEERING] Delle 
RESISTANCE AFTER THE TRANS Dah ieee ie 
WHEN FINDING OPTIMAL PID PARA 
USE JGR seo Oey 
IEF (TIME.LE.1000.) €O TOReeen 


KRESS 
WK ee 
Roe 


k**e ee 


DELSTI= ((XA4*V*V+ XS*R*R+ XO*D* DT ~~ ee eee 


il + DELST1) 
TOTFR2= (X4*V*V + X5*R*R +X1*DU +X2*DU*DU 


+X3*DU*DU*DU + X6*D*D ee 


+ EX - FXX- 63420.) + TOTFR2 


FINDING ACCELERATION VALUES OF U, V;, AND VR “3 


UDOT —=sa ed 
VDOT = (C4*F2-C3*F3)/(C2*e4ees: aan 
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Cke*ks 


C*x*x* 


50 


Ck** 


Cc 


Ce** 


RDOT = (C2*F3-C5*EF2)/(C2*C4-C5*C3) 


WHEN TO PRINTOUT 
IF (ICOUNT.EQ.2) GO TO 50 
COmLems CC 


CONVERT RADIANS TO DEGREES 
YAWDEG= YAW*57.296 


RDEC-R-~5>7.2706 


RDDEG=RDOT*57 .296 
DDEG=D*57 .296 
YAWC=YAWC*57.296 


ONLY USE WHEN GETTING GRAPHS 
MH ( ITIME+1)=DDEG 


T( ITIME+1)=TIME 


ITIME=ITIME+1 
ICOUNT=1 


Essie ir WANT L[O.STOP 


600 IF (TIME.GT.ETIME) GO TO 400 


Se UNTIEGRATION STEP SI4ZE DELT 


DELT=1. 


C INTEGRATION 


ee COULOWINGS ARE FORD TYPE CONTROLLER. 


U=U+QUDOT*DELT 
V=V+VDOT*DELT 
R=R+RDOT*DELT 
YAW=YAW+R*DELT 


Sl=S1+DS1*DELT 


7 


Ree 


kK 


Ree RRA 


KRKEKAKRSK 


S2=S2+(YAWE/KD)*DELT 


FOLLOWINGS ARE FOR PID. 
PIDV2=PIDV2+PIDDV2*DELT 
PIDV3=PIDV3+PIDDV3*DELT 
PIDV4=PIDV4+( YAWE* (PIDK1/PIDT3) ) *DELT 


O° 0 Of: Ooo 


C*** SO MUCH SPEED DECREASING MAKE foE (NONE INEA Ts. 
C*** COEFFICIENT VALUES INCORRECT Gee) rr oars babies ss 


Cees THE SHIP STEBRING erie fig ccs 65 35 
C*** SO RPM (PROPELLER ROTATION PER MINOGE) See 
C*** INCREASED TO PREVENT SHIP SPEEDS Peer. es BSS D0 13 


IF(RPM.GT.95.) GO TO 65 
IF(U.LT.14.7*1.689) RPM=RPM+1./10 


65 IF(U.GT.15.3*1.689) RPM=RPM-1./10. 

@) 

C*** TO FIND SHIP MOTIONS ON THE EARTH SURFACE. ***** 
C*** (REPS. deecoD 1 


XDOT = U*DCOS(YAW) - V*DSIN(YAW) 
YDOT = U*DSIN(YAW) + V*DCOS(YAW) 
X= X 2 por ppEm 
Y= y fo Dol brn 


e 

C***HOW MUCH ADDED FORCES OBTAINED BY INCREASTNG WEP! 
RDIF=4480.*RPM-243908. 
RES=RDIEF-63420. 


C*** WITH USING SHIP ENGINE SPECIFICATIONS, THE FUEL aw 
C*** CONSUMPTION CAN BE FOUND. ee 
EHP=263 . 448*RPM-15153.5312 
IF ( RPM.LE.80. ) SHP = 250*RPM - 11900. 
IF ( RPM.GT.80. ) SHP = 420*RPM - 25500. 


BARFUE = BARFUE + (0.03125*SHP + 95.) * (DELT/86400. ) 
LBFUEL = LBFUEL +(0.000007187*SHP + 0.632625) *SHP 
Cc * (DELT/3600. ) 


Ws 


TIME=TIME+DELT 
ICOUNT=ICOUNT#+1° 
C*** FINDING APPROXIMATE COST FUNCTION wee 
IF(TIME.GT.1000.)ISE=ISE + LAMDA*YAWE**2 
IF(TIME.GT.1000.)ISR=ISR + D**2 
TDIFF=(ISE + ISR) 


é 

GO TO 200 
400 CONTINUE 

©: 

C*** FINDING FUEL CONSUMPTION RATIO or 
DELSTY = DELST1 / (TIME - 1000.) . 
TOTER2 = TOTFR2 / (TIME - 1000.) 

PERC2 = (DELST1/TOTFR2) * 100. 

é 


WRITE(6,1111) TDIFF,DELST1,PERC2 
1111 FORMAT(1X,F7.3,1X,F8.0,1X,F8.3) 


Cc 
C***USE THIS PART TO PRINT ALL COMBINATIONS OF PID***** 
C***PARAMETERS eR KK 
Cc WRITE(8,43)KK1(1),KK2(IK2),TT1(1T1),TT2(1T2) 
Cc Cc 7TTS(1T3),TDIFF 
@ MmGpnem Ss NIC d ) Kie24 1K2),TT1(1T1), TL2(iT2) 
C C /'T3(1T3), TDIFF 
43 FORMAT(' ',5F8.2,2X,F20.5) 
e 
C*** DISSPLA STATEMENTS USING TO GET GRAPHS Aw RK 
XMIN=O 
XMAX=ETIME 


XINC=XMAX/6. 

MAGMAX=-1.E15 

MAGMIN= 1.E15 

DO 150 I=1,ITIME 

IF (MH(1).GI.MAGMAX) MAGMAX=MH(1) 
150 IF (MH(I).LT.MAGMIN) MAGMIN=MH(1) 


ik 


Ss 


C*** USE THIS PART TO FIND OPTIMAL PID PARAMETERS 


€ 
© 
C 
Cos 
C44 
Cee 
C2zZ 
ei 


Gye IIb: 
CALL 
Calg 
Sig: 
CAang 
CAEIE 


Qahe 


CALL 
CALL 
CALL 
CAng 
CALE 
CALL 
CALL 


AREAZD G52 O cme. 
XNAME('TIME (#SEC.&) $',100) 
YNAME('RUDDER (#DEG.&) $',100) 


HEADIN(' So OO ly zaece) 
HEADIN ( ' SEA STATE=8 $',100)1 27a 
HEADIN(' EN. ANGLE=030 SPEED=15 


D TYPE CONT. —SlOO meee s 


HEADIN( ' WF = .53 ( EF =.33) 
S* 1007122 4) 
CROSS 
YINTAX 
GRAF (XMIN, XINC, XMAX, MAGMIN, 'SCALE' , MAGMAX) 


CURVE(T,MH, ITIME,0) 
RESET('CROSS') 
ENDPL(0) 

DONEPL 


KH 


IF (TDIFR ERS Ali eee 


SMALL=TDIFE 


END IF 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 
CONTINUE 

WRITE(8,*) SMALL 


SL Or 
END 
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APPENDIX D 
LIST OF SYMBOLS 


Wo Commanded Heading Angle 


wy Yaw Angle, measured from the vertical xz plane 
to the axes of the ship; positive in the 


positive sense of rotation abiut the z-axes 


We ¥ - We 


Xx Hydrodynamic force components on ship body 


(longitudinal ) 


a6 Hydrodynamic force components on ship body 
(lateral) 
N Resultant total moments acting on a ship 


about the z-axis, yawing moment 


u Velocity components of the origin of the body axes 


relative to the fluid, longitudinal component 


Uy Commanded longitudinal velocity component 
Au eee Uy 

V Transverse component of velocity 

ae Angular velocity of yaw 

u, Vv, r Acceleration components of the origin of the 


the body axes relative to the fluid (longitudinal, 


transverse, and yawing, respectively) 


re) Angular displacement of the rudder, measured from the 


xZ-plane of the ship to the plane of the rudder 


m Mass of ship 


19 


OP 


Coordinate of the center of mass of the ship 


ship relative to body axes 
Values of X, Y, and N at V=r=v=r=0 and u=U, 
Significant wave height 
EXclting force magnaunde 
Exciting force phase angle 
Encounter frequency 
Wave circular frequency 
Encounter angle 
Acceleration of gravity 


Length between perpendiculars 


RHR! 5/2 
Natural frequency of ship's steering 


Distance from the ship pivot point to the orveam 
(QO ..3.0eien) 


Ship length 
Weighting factor 
Amplitude of r, v, and § 

Engine power 

Work done per unit mass (btu / l1bm) 
Fuel mass per unit time (lbm / min) 


wer (ie = 


propeller in calm water case 


t), total ship resistance without 


Thrust deduction ftrvacereon 
Propeller thrust 
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Xo Partial derivative of X with respect to v4 


ae Partial derivative of Y with respect to rv 


Ay Partial derivative of N with respect to v 


XQ YQ.2Q System of reference axes through origin 
of reference axes fixed in the ship 


whose direction fixed in the ship 


wy, Zz System of reference axes whose origin 


and direction remain fixed in the earth 
Notes: 


[ee oem cor cl oucections, forces, distances, veloci- 
ties, and accelerations are positive downward along the 
Z-axes, positive to starboard along the y-axes, and positive 


forward along the x-axes and similarly along the Xj Yo and 
“O0-axes 


2. Signs of all angles, angular velocities, angular 
accelerations, and moments are positive if clockwise when 


facing in the positive direction of appropriate axes. 
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